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We study the electric quadrupole (E2) contribution to resonant x-ray scattering (RXS). Under the 
assumption that the rotational invariance is preserved in the Hamiltonian describing the intermediate 
state of scattering, we derive a useful expression for the RXS amplitude. One of the advantages 
the derived expression possesses is the full information of the energy dependence, lacking in all the 
previous studies using the fast collision approximation. The expression is also helpful to classify 
the spectra into multipole order parameters which are brought about. The expression is suitable 
to investigate the RXS spectra in the localized / electron systems. We demonstrate the usefulness 
of the formula by calculating the RXS spectra at the Ce £2,3 edges in Cei-^La^Be on the basis of 
the formula. We obtain the spectra as a function of energy in agreement with the experiment of 
Ceo.7Lao.3B6. Analyzing the azimuthal angle dependence, we find the sixfold symmetry in the a-a' 
channel and the threefold one in the ct-tt' channel not only in the antiferrooctupole (AFO) ordering 
phase but also in the antiferroquadrupole (AFQ) ordering phase, which behavior depends strongly 
on the domain distribution. The sixfold symmetry in the AFQ phase arises from the simultaneously 
induced hexadecapole order. Although the AFO order is plausible for phase IV in Cei-^La^Bg, the 
possibility of the AFQ order may not be ruled out on the basis of azimuthal angle dependence alone. 

PACS numbers: 78.70.Ck, 75.25.+Z, 75.10.-b, 78.20.Bh 
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I. INTRODUCTION 

Resonant x-ray scattering (RXS) has recently at- 
tracted much attention, since strong x-ray intensities 
have become available from the synchrotron radiation. 
It is described by a second-order process that a core elec- 
tron is excited into unoccupied states by absorbing inci- 
dent x-rays and that electron is recombined with the core 
hole by emitting x-rays. The RXS has been recognized as 
a useful probe to investigate spatially varying multipole 
orderings, which the conventional neutron scattering is 
usually difficult to detect. 

For probing the spacial variation of order parameters, 
x-ray wavelengths need to be order of the variation pe- 
riod. In transition metals, the K edges in the dipole 
(El) transition are just fitting for this purpose. Actu- 
ally, by using the K edge, the possibility of the orbital 
ordering has already been explored in transition-metal 
compounds^ 2 - The RXS intensities are observed at su- 
perlattice Bragg spots, which are interpreted as originat- 
ing from the modulation in the Ap band, since the process 
involves the excitation of a Is electron to unoccupied Ap 
states. 

Because the ordering pattern is usually controlled by 
electrons in the 3d band, the mechanism which causes the 
modulation is not necessarily trivial. Actually, for most 
of transition- metal compounds, both experimental stud- 
ies and theoretical studies based on electronic structure 
calculations have revealed that the RXS intensities are 
brought about by the hybridization between the Ap band 
and the 2p band of the neighboring anions rather than 



the direct Coulomb interaction between the electron in 
the Ap band and electrons in the 3d band. 3,4 This result 
is reasonable because of the extended nature of the Ap 
state. 

On rare earth metal compounds such as CeBg, 
DyB2C2, the £2,3 edges in the El transition are used 
because of the requirement for x-ray wavelength ^SiS 
The RXS spectra in the El transition from the antifer- 
roquadrupole (AFQ) phase of CeB 6 were studied both 
experimentally^ and theoreticallyiSiiS Although the ex- 
periments and the theory give sufficiently consistent re- 
sults, the relation to the multipole orderings which 4/ 
electrons mainly involve is rather indirect, since the res- 
onance is caused by the excitation of a 2p electron to 
5d states. This shortcoming may be overcome by using 
the quadrupole (E2) transition at the £2,3 edges, where 
a 2p electron is promoted to partially filled 4/ states. 
Using the E2 transition has another merit that octupolc 
and hexadecapole orderings are directly detectable. This 
contrasts with the El transition, where only dipole and 
quadrupole orderings are detectable. Of course, intensi- 
ties in the E2 transition are usually much smaller than 
those in the El transition. 

In this paper, we derive a general formula of the RXS 
amplitudes in the localized electron picture, in which the 
electronic structure at each atom is assumed to be well 
described by an atomic wavefunction under the crystal 
electric field (CEF). Historically, the research in such 
a direction was started by using the framework bor- 
rowed from resonant 7-ray scattering^! Starting from 
the works by Blume and Gibbs 12 and by Hannon et ali&, 
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the form of the RXS amplitude had been investigated in 
several works ^i^igiil The RXS amplitude can be sum- 
marized into an elegant form by using vector spherical 
harmonics. Unfortunately, it has little practical usage 
because it is difficult to deduce meaningful information 
when there is no restriction on the intermediate state 
of the scattering process. A widely-adopted approxima- 
tion for practical use is the so-called "fast collision (FC) 
approximation" . This replaces the intermediate state en- 
ergy in the energy denominator of the RXS amplitude by 
an averaged value, allowing the denominator out of the 
summationiiiiiiSiiSiii Thereby, the multiplet splitting 
of the intermediate state is neglected, leading to an as- 
sumed form (usually a Lorentzian form) for the energy 
profile. 

However, recent experiments show deviation from the 
Lorentzian form in several materials pi^iiS We improve the 
situation by taking the energy dependence of the interme- 
diate state under the assumption that the intermediate 
Hamiltonian describing the scattering process preserves 
spherical symmetry. This assumption is justified when 
the CEF energy and the intersite interaction are much 
smaller than the multiplet energy in the intermediate 
state as is expected in many localized electron materi- 
als. We have already reported the formula for the El 
transition, having successfully applied to the analysis of 
the El RXS spectra in URu 2 Si 2 and NpO 2 , 20 i 21 i 22 This 
paper is an extension of those works to the E2 transi- 
tion. The obtained formula makes it possible to analyze 
the energy profiles of the spectra in contrast with the FC 
approximation. In addition, the formula is suitable to 
discuss the relation of the RXS spectra to multipole or- 
der parameterSfiSiiSiSSiSi 2 ^ 2 ^ because it is expressed by 
means of the expectation values of the multipole order 
parameters. 

We demonstrate the usefulness of the formula by cal- 
culating the RXS spectra in multipole ordering phases of 
Cei_ x La x Bg. First, we investigate the E2 RXS spectra 
at the Ce L 2 ,3 edges from the AFQ ordering phase (phase 
II) in the non-diluted material CeBV Analysis utilizing 
our formula reveals that the E2 RXS spectra in phase 
II consist of a mixture of the quadrupole and hexade- 
capole energy profiles. The calculated intensities suggest 
the possibility that he E2 signal at the Ce L 2i 3 edges can 
be detectable in this material. 

For the intermediate doping range x rs 0.3 ~ 0.5, 
Cei-^La^Bg falls into a new phase (phase IV) whose 
primary order parameter is not well established yet^i 
Recently, the E2 RXS signals at the Ce L 2 edge have 
been detected for an x = 0.3 sampled From the az- 
imuthal angle (ip) dependence of the peak intensity, it 
was claimed that the antiferrooctupole (AFO) ordering 
phase is the most probable candidate because the sym- 
metry of the i/'-dependence, sixfold and threefold in the 
a-a' and a-n' channels respectively, is deduced from the 
theory in good agreement with the experiment. 25,28 How- 
ever, the relative intensity between two channels depends 
strongly on the domain distribution, and deviates about 



factor two from the experimental one if the contribution 
from four domains are added with equal weight. The 
origin of this discrepancy is still unanswered. It will be 
pointed out that the RXS peak intensity from the AFQ 
phase concomitant with the induced hexadecapole con- 
tribution also gives rise to the same symmetry of the 
^-dependence as obtained from the AFO phase. Thus, 
although the AFO order is plausible in many respects, it 
seems difficult to rule out the AFQ order on the basis of 
the azimuthal angle dependence alone. In addition, we 
calculate the energy dependence of the RXS spectra at 
the Ce L 2> 3 edges. Assuming both the AFO and AFQ or- 
ders, we obtain the spectral shapes at the L 2 edge, which 
agree with the experimental one for Ceo.rLao.aBgiS 5 . On 
the other hand, the spectral shapes at the L3 edge are 
found slightly different between two phases, with inten- 
sities the same order of magnitude of the reported one at 
the L 2 edge. 

The present paper is organized as follows. A general 
expression for RXS amplitudes is obtained in Sec. ITU 
Analysis of RXS spectra in Cei-^La^Bg are presented in 
Sec. IIIII Section llVl is devoted to concluding remarks. In 
Appendix, we show several expressions required to obtain 
the RXS amplitude formula. 



II. THEORETICAL FRAMEWORK OF RXS 



A. a second-order optical process 



The RXS is described by a second-order optical pro- 
cess, where a core electron is excited to unoccupied states 
by absorbing x-rays and that electron is recombined with 
the core hole by emitting x-rays. Since the wavefunction 
of core electron is well localized, the RXS amplitude may 
be given by a sum of contributions from individual ions. 
Using a geometrical arrangement shown in Fig. ^ we 
express the RXS amplitude /(e, e', k, k', w) for the inci- 
dent x-ray with momentum k, polarization e, and the 
scattered x-ray with momentum k', polarization e' as 

/(e,e',k,k» = ^/W(e,e',k,k», (2.1) 

/W(e,e',k,k» oc _L]Te- 4G ^ 

j 

x MW(e,e',k,k», (2.2) 



3 



a 

scattered photon 




sample 



FIG. 1: Geometry of the RXS experiment. Photon with 
polarization a or ir is scattered into the state of polarization 
a' or 7r' at the Bragg angle 9. The azimuthal angle ip describes 
the rotation of the sample around the scattering vector G. 
The (112) surface is in the scattering plane at = 0. 
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hw - (E A - E ) + iT 
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Mj 2) ( e , e ',k,k',c) - T ^^(k',e')v(k,e 

(^o|^|A)(A|v,j#o) 
Huj-(E a - E ) + iT 
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where G (= k' — k) is the scattering vector, and N is 
the number of sites j's. The jV'o) represents the ground 
state with energy Eq, while |A) represents the inter- 
mediate state with energy E\. The T describes the 
life-time broadening width of the core hole. Equation 
(|2.3() describes the El transition, where the dipole op- 
erators a^j's are defined as X\,j — Xj, X2.j = Vj, and 
j = Zj in the coordinate frame fixed to the crys- 
tal axes with the origin located at the center of site j. 
Equation (|2.4() describes the E2 transition, where the 

quadrupole operators are defined by Zij — -^{x" 2 - " 2 
%2,j = §(3z| - Tj), h,j = V^yjZj, z 4J = \/3. 
Z5,j — V^XjUj. Factors g p (k, e) and g M (k',e') with 

k = k/|k| and k = k'/|k'| are defined as a second-rank 
tensor, 



ZjXj and 



^(A,B) 



& {A X B X - A y B y ) for/i=l, 

\ {3A Z B Z - A • B) for /i = 2, 

^(A y B z + A z B y ) for /i = 3, 

+ for /i = 4, 

& (A x B y + A y B x ) for /i = 5. 



(2.5) 



Note that the quadrupole operators are expressed as 
z^-j — Qfj.( r ji r j)- 



B. Energy profiles 

In localized electron systems, the ground state and the 
intermediate state are described in terms of the eigen- 
functions of the angular momentum operator, | J, to), at 
each site. At the ground state, the CEF and the intersite 
interaction usually lift the degeneracy with respect to to. 
We write the ground state at site j as 



]cj(m)\J,m) 



(2.6) 



In the intermediate state, however, the CEF and the in- 
tersite interaction may be neglected in a good approxima- 
tion, since their magnitudes of energy are much smaller 
than those of the intra-atomic Coulomb interaction and 
the spin-orbit interaction (SOI) which give rise to the 
multiplet structure. Thus the Hamiltonian describing 
the intermediate state is approximated as preserving the 
spherical symmetry. In such a circumstance, the inter- 
mediate states are characterized by the total angular mo- 
mentum at the core-hole site, that is, |A) = \J',M,i) 
with the magnitude J' and the magnetic quantum num- 
ber M. The corresponding energy is denoted by Ejij, 
where we introduce the index i in order to distinguish 
multiplets having the same J' value but having different 
energy. 

In the following, we discuss only on the E2 transition 
(Eq. (|2.4[l ). because the El transition has been fully ana- 
lyzed in our previous paper. 21 First, we rewrite Eq. I|2.4|) 
as 

Aff( e ,e',k,k',c) = ^^(k%>^k,e)M^(^), 



ii, fi' 



J'M.i 



x (J',M, i|v,j#o), 



(2.7) 



(2.8) 



with 



EAw,J') 
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hw - (Ej, ti - E Q )+iV 



(2.9) 



Then, inserting Eq. (|2.tj|) for the ground state into 
Eq. (|2.8fl . we obtain 



M (2) 



>(?» = £ c*(m) Cj ( m ')M^r m M, (2.10) 



with 



J' 1=1 



J' 



J2 (J,m\^\J',M,%)(J , ) M,i\z ft ,\J,m'X2.U) 



M=-J> 
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We have suppressed the index j specifying the core-hole 
site. The number of the multiplets having the value J' is 
denoted by Nji . The selection rule for the E2 transition 
confines the range of the summation over J' to J' = 
J, J±l, J ±2. 

Now we analyze the matrix element of the type 
(J, m\Zfj,\J' 7 M) by utilizing the Wigner-Eckart (WE) the- 
orem for a tensor operator^ 



(J,m\ Stt \J'M) 



J' +1X1-1 



-1) 

J' 2 
M n 



V2J+ 1 



J 



(J||F 2 ||J'),(2.12) 



with s±2 — (z\ ± iz5)/V2, s±i — +(24 ± i^3)/V2 and 
s o = z 2 - The symbol (J||V2||J') denotes the reduced 
matrix element of the set of irreducible tensor operator of 
the second rank. Because of the nature of the quadrupole 
operators, a condition \m — ml\ < 4 has to be satisfied for 

non- vanishing M^J n,m (u>). After a straightforward but 
tedious calculation with the help of the WE theorem, we 
obtain non-zero M^ m,m (w)'s. Then, we perform the 
summation over m and m! in Eq. (|2.1U|) . The result is 
summarized by introducing the expectation values of the 
components of the multipole operators as follows: 



2v+l 



M,M'(V>0|4 |V>0 



(2.13) 

where the A th component of rank v tensor zf in real 
basis (1 < A < 2v + 1) is defined in Table [Q The 

(n) 

is constructed from the irreducible spherical tensor 
through the unitary transformation U^ u > . The definitions 
of Tn^ and as well as the energy profile a^\{u)) are 
given in Appendix ^ The matrix element of is 
expressed as 



(2||r„||2) V 5 



TABLE I: Definition of the operator equivalence of the multi- 
pole order components. The overline denotes the symmetriza- 
tion, for instance, ~XW = X 2 Y + XYX + YX 2 . 

z< 1} = J, 

= j. 

— IVT 



f =O az2 _ r 2 = ±[3J 2 -J(J +1)] 
3 ^ys ~2~\JyJz ~t" JzJy\ 



= O 



_ n/3 



J x J y ~~\~ J 1) J x 



'■ T S 
■■ T° 



s/15- 



-1 

>> 

8 (3) =T a = 1 [2J 3 

,(3) *Tnj3 VT5 

'6 " J » " 6 

.(3) 



J x J y J z 

|[2Jg- ^(Jg_+Jg) ] 

■M-g + Jg)] 



^ 2 ) 



jg) 



(4) 



r (4) 



= HJ = 



; i [j; +Jt+Jt -§./(./ +i) {./(./ 1 }j 

^ [So/i-W - { J(J + 1) - 1} (4 - Jy)} 

< 4 > = Hi = $ [35Ji - 30.7(7 + 1)J 2 + 3J( J + 1)[J(J + 1) - 2] 
+25 - 7[4 + J 4 - JpJ]] 



,(4) 
s 4 

(4) 
H 
.(4) 



4 4) = fff = ^[2J*J V J. ~ ( Jg J» + J»J|) ] 

(JJJ. + J 2 J|)] 



= H.' 3 = 



Zq — — o J? J x J u 



Eq. 12.14|l . we obtain the final expression, 



fc 2 4 



2i/+l 



9^ a BaH U Pr(^e',k ) k)(Vol4 !/) lV'o), 

i^=0 A=l 



£=-2l' = -2 



(2.16) 



with 



x E ^WV)*[t/ (2)t W' 

n— — v 

= (-TV2I7TT E E {-Yu$ 

t=-2i' = -2 



(2||T,||2) 



(5 + i/)! 



2^ V 5(4- j/)!' 



(2.14) 



(2.15) 



where -Pj^'s are the geometrical factors defined as 



= v^ttt e (-)"^ E 

n=—v 

x (? m (e',k )q„_ m (e,k 



m— — u 



2 2 v 

m n — m —n 



(2.17) 



Finally, inserting Eq. I|2.13|) into Eq. I|2.7|) and using 



Those for v — 0, 1 and 2 are expressed as relatively 
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simple forms: 
P 1 (0) (e ) e',k,k / ) 

FW(e,e',k,k) 



V5 



(k ■k)(6'-6) 



1 



(e' ■ e)(k x k) A 



+ (k'-e)(e'xk) M + (e'-k)(k x e) M) 
Pf (6,6',k,k) = -^[(e'-e^k) 

+(k ■ k)^(e, e') + g M (k x k, e' x e) 



+ (k -6)(e'.k) , 
(2.18) 

(k-k)(e' x e) M 

(2.19) 



(2.20) 



For v = 1, indices fx = 1, 2 and 3 serve as the Cartesian 
components and z, respectively. The corresponding 
expression of Pjf^s for v = 3,4 have complicated forms, 
which are summarized in Appendix B. 

An expression similar to Eq. I|2.16[) has been derived by 
the FC approximationiiiii!ii£ii£ii£ However, this scheme 
has to put by hand the energy dependence. The present 
theory gives an explicit expression of the energy depen- 
dence, which is separated from the factor relating to the 
order parameter. Thus, the choice of the CEF param- 
eters in the ground state does not affect the shape of 
energy profiles a^\{ui). 



III. APPLICATION TO MULTIPOLE 
ORDERING PHASES IN Cei-^La^Be 

In this section, we demonstrate the usefulness of 
Ea. (|2.16|l by analyzing the RXS spectra in the E2 tran- 
sition at the Ce £2,3 edges from Cei-^La^Bg. 



A. Phase II in CeB 

The parent material CeBg experiences two-step phase 
transitions. It undergoes the first transition from param- 
agnetic (phase I) to an AFQ state (phase II) at Tq = 3.2 
K and the second transition to an antiferromagnetic 
(AFM) state (phase III) at T N = 2.4 K under no ex- 
ternal magnetic field. The AFQ order is known to be a 
Neel-type with a propagating vector Q = 

These phase transitions have been theoretically stud- 
ied in a localized electron scheme, where each Ce ion is 
assumed to be trivalent in the 4/ 1 configuration. Its 
ground multiplet is a Tg quartet confined within the 
J = § subspace under the cubic symmetry. Using states 



{\J Z = to)}, the four bases |±,er) (cr 
pressed as 

_ 5 

2 




I- T> = 



, l) may be ex- 

(3.1) 
(3.2) 



and |±,|) by replacing \m) with | — to). The intersite 
interaction may lift the fourfold degeneracy, leading to 
multipole orderings. Shiina et al. have derived such inter- 
action from a microscopic model and obtained the phase 
diagram in agreement with experiments^2*2i 

Instead of pursuing this direction, we simply assume 
the ordering pattern, and calculate the RXS spectra. 
The assumed ordering pattern selects a particular en- 
ergy profile according to Eq. (|2.16() . Note that the quar- 
tet Tg consists of 16 degrees of freedom, which are ex- 
hausted by three components of dipole, five components 
of quadrupole and seven components of octupole opera- 
tors as well as an identical operator. Thereby the hcxadc- 
capole operators , i?| , H\ , H$ , H@ and are equiva- 
lent to identical operator, O x 2_ y 2 1 3z 2_ r 2, O yz , O zx and 
O xy , respectively, while H^ y z = 0. Therefore, as long 

as a contribution from a~^\(ui) exists, that from a^\(uj) 
automatically exists. 

The order parameter in phase II is believed to be the 
C^y-type. Operator O xy has two degenerate eigenstates 
of eigenvalue — 1 and two degenerate eigenstates of eigen- 
value +1, that is, 



-10 
0-100 
10 
1 



(3.3) 



within the bases of eigenfunctions. The AFQ phase may 
be constructed by assigning two degenerate eigenstates 
with eigenvalue —1 to one sublattice and those with 
eigenvalue +1 to the other sublattice. The degeneracy of 
the Kramers doublet would be lifted in the AFM phase 
with further reducing temperatures. Within the same 
bases in order, typical dipole and octupole operators are 
represented by 



■h. = 



( 



9 -I \ 



7 

3 6 





1 1 



i3V5 \ 





-i3\/5 

tiVE 

V -i3\/5 / 



(3.4) 



(3.5) 



These forms indicate that the O xy order could accompany 
neither the J z order nor the T£ order. Therefore, the 

O xy order selects the energy profiles a^\{uj) and a^\{uj) 
according to Eq. (|2.1t)|) . 

In the actual calculation of energy profile, we take into 
account full Coulomb interactions between 2p and 4/ 
electrons, between 2p electrons, and between 4/ electrons 
in the configuration (2p) 5 (4/) 2 . The spin-orbit interac- 
tion (SOI) of 2p and 4/ electrons are considered too. 
The Slater integrals necessary for the Coulomb interac- 
tions and the SOI parameters are evaluated within the 
Hartree-Fock approximation. 32-33 
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FIG. 2: The RXS spectra at the Ce L 3 (top) and L 2 (bottom) 
absorption edges from the AFQ phase (phase II). F — 2.0 eV 
and 1.0 eV, G = (§§§), and ip — 0. Only the spectra in the 
a-a' channel are displayed. 



Figurc[21shows the RXS spectra as a function of photon 
energy, calculated with the core-hole lifetime width T = 
2.0 eV and 1.0 eV. The energy of the Ce 2p-core level is 
chosen such that the peak of the RXS spectra at the Ce 
L3 edge coincides with the experiment for CeBg. We find 

that the absolute value of a^\(uj) is much smaller than 

(2) 

that of a E 2{uj). However, the smallness is compensated 

by a large value of (ipo\z^ \ipo) , and thereby both terms 
contribute to the intensity. The calculated spectra show 
asymmetry and some structures, which depend on the T 
value. It may be appropriate to use T = 2.0 eVi^i 

When the O xy order is realized, the O yz and O zx or- 
ders are also possible to be realized. In actual crystals, 
three orders may constitute domains, whose structure af- 
fects the azimuthal angle dependence of the RXS spectra. 
Figure 13 shows the peak intensity as a function of %j) for 
the scattering vector G = (§§§)■ The origin of tp is de- 
fined such that the scattering plane includes the a-axis. 
It depends strongly on domains. An incoherent addition 
over the contributions from three domains is performed. 
In the a-o' channel, the term of v — 2 in Eq. I|2.16|l is in- 
dependent of ij) so that the sixfold symmetry comes from 
the term of v = 4. On the other hand, the threefold 
symmetry in the cr-7r' channel arises from both the term 
of v = 2 and that of v = 4. 



B. Phase IV in Cei-^La^Be 

The La diluted material Cei-^La^Bg with x ~ 0.3 ~ 
0.5 has an additional phase IV whose order parameter is 
not well understood yet^i Although a large discontinu- 
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FIG. 3: The peak intensities of RXS as functions of az- 
imuthal angle at the Ce L3 edge from the AFQ pahse (phase 
II). r = 2.0 eV. Three domains are assumed to have equal 
volumes. The solid and broken lines represent the a-a' and 
<j-7r' channels, respectively. 



ity in the specific heat curve suggests the existence of the 
long range order no neutron scattering experiment has 
found an evidence of long range magnetic order It 
is suggeste d 38 ! 39 that the AFO order characterizes phase 
IV, which is supported by the observation of the trigo- 
nal distortion. 40 Recently, Mannix et al. measured the 
RXS spectra at the L2 edge in the E2 transition in 
Ceo.7Lao.3B6, claiming that the signal arises from the 
AFO orderiSi The analysis of the azimuthal angle de- 
pendence by Kusunose and Kuramoto supports the AFO 
order in phase IV. 28 However, there exists at least one 
prominent discrepancy between the experiment and the 
theory about the azimuthal angle dependence which we 
shall address later. 

Keeping two possibilities, the quadrupole and octupole 
orders, for phase IV, we analyze the spectra on the ba- 
sis of Eq. (|2.16|l . Since the trigonal distortion is observed 
along the body-diagonal direction, we assume that the or- 
der parameter is of lf u type (lf u = {T^+T^+T^)/V3) 

or O m type (O m = (O xy + O yz + O zx )/V3). The T a 
type can be ruled out because this type carries a sub- 
stantial antifcrromagnetic moment, which is against the 
experimental finding. Since [Ty^Om] = 0, both oper- 
ators are simultaneously diagonalized. Within the bases 
of eigenfunctions, they are represented as 



111 



Ox 
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(3.6) 



(3.7) 
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Within the same bases, the dipole operator J m (= (J x + 
J y + Jz)l v3) is represented as 

/ zx \ 

j _ *i* o 

Jm- _7 , (3.8) 

V o o o I) 

where z x = *§(1 + ill-y/2). 

_Z. AFO order 

The AFO order may be constructed by assigning the 
eigenstate of the Tf n with eigenvalue — 3\/T0 to one 
sublattice and that with 3VT0 to the other sublattice. 
Then, the order parameter vector ({Tg), {T£), (Tg )) is 
pointing to the (111) direction. Equation l|3.8[l indi- 
cates that the AFO order could carry no magnetic mo- 
ment, which is consistent with the experiment. Equa- 
tion H3.7fl indicates that the AFO order accompanies the 
ferroquadrupole order, not the AFQ order. Therefore, 
according to Eq. i|2.16[) . the RXS energy dependence is 

purely characterized by |a£2(w)| 2 . 

Figure |1| shows the calculated la^^w)! 2 as functions 
of the incident photon energy u) at the Ce L2 and L3 ab- 
sorption edges with T = 2.0 eV and 1.0 eV, in comparison 
with the experiment of Mannix et al (the non-resonant 
contribution is subtracted from the data)^ In the cal- 
culation, we use the same Slater integrals and the SOI 
parameters as in phase II. The spectral shapes depend 
strongly on the absorption edge they are observed. In 
particular, the tail part of the spectra at the L3 edge is 
drastically different from that at the L2 edge. This fact 
might be helpful to identify the character of the ordering 
pattern if the spectrum at the L3 edge is experimentally 
available. Since the peak intensity at the L 3 edge is about 
20 % of that at the L2 edge, it can be said that experi- 
mental observation has a legitimate chance at the former 
edge. The L2 spectral shape reproduces well the exper- 
imental one showing broad single peak structure with a 
hump in the high energy region. 

The energy profile la^^o;)! 2 looks similar to the spec- 
tral shape in phase II (Fig. [5J for T = 2 eV. One differ- 
ence is a dip found at the L 3 edge in la^Cw)] 2 , which 
is absent in Fig. If the T is as small as 1 eV, the 
differences are emphasized around the tail part of the 
high energy region, because multiplet structures of the 
intermediate state are emphasized, 20 Note that, although 
|a£2((xj)| 2 is about two order of magnitude smaller than 

|oj^2( w )| 2 i the smallness is compensated by the large fac- 
tor of \{T,P y _ z )\ 2 ~ 90, resulting in the same order of 
magnitude of the spectral intensity as in phase II (Fig. 
0). 

If the octupole order parameter vector 
((Tf), (T^>, (If)) can point to the (111) direction, 
it is also possible to point to the (111), (111), and (111) 




5705 5715 5725 5735 5745 




6140 6160 6180 

a) [eV] 



FIG. 4: Energy profile |a^(aj)| 2 at the Ce L 3 (top) and 
L2 (bottom) absorption edges at G = (§§§)■ Curves with 
r = 1.0 eV (broken lines) are multiplied by factors 0.29 at 
the L3 edge and ... at the L2 edge to have the same peak 
intensities as those with F = 2.0 eV (solid lines). Open circles 
are experimental data in Ceo.7Lao.3B6, in which non-resonant 
contributions are subtracted as explained by the authors.—- 



directions. These four orders usually constitute domains. 
The azimuthal angle dependence is different for different 
domains, as shown in Figs. EJa) and (b). If you collect 
the contributions from domains with equal weight, the 
maximum intensity in the a-ir' channel become nearly 
equal to that in the a-a' channel. The experimental data 
show that the maximum intensity in the a-ir' channel is 
about the half of that in the a-a' channel, as shown in 
Fig.[5Jc). This may be attributed to the slightly different 
setup for different polarizations and/or to the extrinsic 
background from the non-resonant contribution, as 
discussed by Kusunose and Kuramoto— - They reduced 
the intensity in the a-ir' channel by simply multiplying 
a factor 0.6. Another possibility is that domain volumes 
are different among four domains. Collecting up the 
contributions with ratio 3:1:1:1 from the (111), 
(111), (111), and (111) domains, we have the result 
similar to that simply multiplying a factor 0.6 to the 
intensity in the a-ir' channel, as shown in Fig. [5fc). 
Thus, the sixfold and threefold symmetries in the a-a' 
and <7-7r' channels are well reproduced in comparison 
with the experiment. 



2. AFQ order 

The AFQ order may be constructed by assigning the 
eigenstates of Om with eigenvalue —1 to one sublattice 
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FIG. 5: Peak intensities of the RXS spectra at the Ce L2 
edges in the AFO phase as functions of azimuthal angle. Pan- 
els (a) and (b) display the peak intensities in the a-a' and the 
a-iv' channels, respectively, where the solid (Dl), broken (D2), 
dotted (D3), and broken-dotted (D4) lines represent the peak 
intensity of the domains (111), (HI), (HI), and (111), re- 
spectively. Panel (c) shows the intensities collecting the con- 
tributions from the domains (111), (HI), (111), and (111) 
with ratio 3:1:1:1. The solid and broken lines represent the 
intensities in the a-a' and a-n' channels, respectively. Filled 
and open circles are the experimental data for Ceo.7Lao.3B6P— 



FIG. 6: Peak intensities at the Ce L2 edges in the AFQ phase, 
as functions of azimuthal angle. Panels (a) and (b) display 
the peak intensities in the a-a' and the a-n' channels, respec- 
tively, where the solid (Dl), broken (D2), dotted (D3), and 
broken-dotted (D4) lines represent the peak intensity of the 
domains (111), (HI), (HI), and (111), respectively. Panel (c) 
shows the result collecting the contributions from the domains 
with equal weight. The solid and broken lines represent the 
a-a' and a-n' channels, respectively. The curve in the latter 
channel is multiplied by a factor 0.6. Filled and open circles 
are the experimental data for Ceo.7Lao.3B6i— 



and those with +1 to the other sublattice. The AFQ 
order accompanies no AFO order. The difference from 
phase II is that the order parameter ((O yz ), (O zx ), (O xy )) 
is pointing to the (111) direction. Therefore, the spec- 
tral shape as a function of energy is nearly the same as 
in phase II. The azimuthal angle dependence depends 
strongly on domains, which is shown in Figs. Ufa) and 
(b). The sixfold symmetry in the a-a' channel mainly 
comes not from the (111) domain but from the other do- 
mains, since the contribution from the former domain is 
constant with varying the azimuthal angle. Collecting 
the contributions from four domains with equal weight, 
and reducing the intensity in the a-n' channel by mul- 
tiplying a factor 0.6 in the same way as Kusunose and 
Kuramoto adopted^ we obtain the result in agreement 
with the experiment at least in a symmetrical point of 
view. Although the amplitude of the oscillation in the a- 
a' channel is too small compared with the experimental 
one, the situation may be changed if the subtraction of 
the non-resonant contribution in the a-a' channel and/or 
that of the enigmatic El contribution in the a-n' channel 
are/is reevaluated. Actually, the discrepancy about the 



relative intensity between two channel may be attributed 
to the consequence of these subtraction process. 

We now turn to our attention to the energy dependence 
of the spectra. Owing to our formula Eq. (|2.1tj|) . the 
spectral shapes from the AFQ phase with Om type in 
phase IV are the same as those with O xy type in phase 
II (Fig. |2J). Therefore, the energy dependence at the L2 
edge is similar to that obtained from the AFO phase. 
On the other hand, the energy dependence at the Lsedge 
from the AFQ phase is different from that obtained from 
the AFO phase, which may help the identification of the 
ordering pattern realized in this material. 



IV. CONCLUDING REMARKS 

We have derived a general formula of the RXS ampli- 
tude in the E2 transition. The derivation is based on the 
assumption that the Hamiltonian describing the interme- 
diate state of the scattering process preserves the spher- 
ical symmetry. The obtained formula is applicable to 
many / electron systems where a localized scheme gives 
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a good description. Although similar formulae have al- 
ready been obtained^*^^*^ 4 ^ the present formula has 
two prominent advantages. One is that it is able to cal- 
culate the energy profile of the RXS spectra, because our 
treatment is free from the fast collision approximation 
adopted in the previous works. The other is that it is 
conveniently applicable to the systems possessing multi- 
pole order parameters. 

We have demonstrated the usefulness of the de- 
rived formula by calculating the E2 RXS spectra in 
Cei_ x La x B6. Phase II is believed to be an AFQ or- 
der of O xy type, and our formula dictates that the en- 

(2) 

ergy dependence is given by a combination of a E 2(u>) and 

oS E2 {oj). We have obtained the RXS intensity in the same 
order of intensity as obtained by assuming the AFO or- 
der. This suggests that the E2 signal is detectable from 
phase II, although only the El signal has been reported 
in phase II of CeB 6 . 5,6 Subsequently, we have calculated 
the RXS spectra by assuming the T-f 1:L -type AFO order, 
in order to clarify the order parameter of phase IV. The 

energy dependence \a E l ( w ) 1 2 nas been obtained at the L 2 
edge in agreement with the experiment in Ceo.7Lao.3Bg. 25 
Unfortunately this is not used to discriminate between 
the AFO and AFQ orders, because the spectral shapes 
are nearly the same in the two ordering phases. On the 
other hand, the spectral shape at the L3 edge has been 
found slightly different from the L2 edge, which might 
help the identification of the ordering pattern. For the 
azimuthal angle dependence, we have reproduced the six- 
fold and threefold symmetries by assuming the AFO or- 
der, in agreement with the previous theoretical analy- 
sis and the experiment The intensity in the a-w' 
channel becomes nearly equal to that in the a-a 1 chan- 
nel with the equal volume for four domains, while in the 
experiment the intensity in the former channel is found 
nearly half of that in the latter. This discrepancy may 
be removed by assuming uneven volumes among four do- 
mains. We have also analyzed the azimuthal angle de- 
pendence by assuming the Om-type AFQ order. It is 
found that the simultaneously induced hexadecapole or- 
der gives rise to the sixfold and threefold symmetries. 
Although the agreement with the experiment is quan- 
titatively not good, it may be difficult to rule out the 
AFQ order from phase IV on the basis of the azimuthal 
angle dependence alone. Since it depends strongly on the 
domain distribution, experiments controlling the domain 
distribution, if possible, might be useful to clarify the 
situation. 
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TABLE II: Irreducible tensor operator T„ with the spherical 
basis. 



rank 






V 


n 




1 


± 1 






2 


± 2 


1 / 3 t2 

2 V 2 J ± 




± 1 









\ [3 J! -./(./ + 1)] 


3 


± 3 
± 2 
± 1 


=F^j J±[15J 2 2 ± 15 J z - 3J(,7 + 1) + 6] 







i[5J^ - 3J(,/+ l)J Z + J Z ] 



±3 TtFJ±(2./z ±3) 

±2 [7 J 2 ± 14J Z - J(J + 1) + 9] 

±1 T^r J ±l UJ l ± 21J l + 19J z 

-6J(J + 1) J z =p 3,7 (J + 1) I 6] 
I [35 J* - 30 J(J + 1)J| + 25 

+3J 2 (./ + 1) 2 -6J(J + 1)] 



APPENDIX A: DEFINITIONS OF SOME 
QUANTITIES USED IN SEC. El 

Let us define irreducible tensor operator of rank v with 
the spherical basis. The n-th component (— v < n < v) 

Tn^ is defined recurrently as 

^ - <-»"/W J - (A1) 

[j_,rM] = + + (A2) 



Expressions for Tn s are listed in Table ITT1 up to rank 
four. We can find (2v + 1) x (2v + 1) unitary matrix 
which connects the tensor operator with the spherical 
component T« and that with the Cartesian component 
which satisfies 



*i" ) = £ E & ,2 f ) ' (A3) 

n— v 

and inversely, 

TM= EV^U^- (A4) 

A=l 

Explicit form of U^ v > is summarized in Table ITTT1 

Finally, we show the explicit forms of the functions 
a'^\(uj), which give the energy profiles coupled to the 



10 



TABLE III: Unitary matrix which connects the tensor op- 
erator with the Cartesian basis and that with the spherical 
basis. 

l 














I 75 




75 





73 







1 





V (2) 



-.1 







V V2 






7f_ 

73 
o 






-.2 





,7(3) 





_ v2 

4 

-4i 



4 





V30 
12 


12 



~75 
o 
o 
o 



V2 










1 





vl 

4 

_ v3; 
4 



V5 
4 

4i 





{) 



1 {) 
1 
{) 






o 

_ vS 

4 



_ vl 

4 





7^\ 





75 / 





4 










4 





-s 



V2 




{) 





"71 



V21 
6 



•STE 
6 











v 2 




{) 





1 

_2_ 



-1 



\ 

yl 
i 

_ vS; 

4 



_ V3 

4 









i 

(I 



-7= 



V2 









75 





_ v5 

4 





%/30 \ 
12 \ 



y/42 
12 




75 






The energy dependence is contained in the functions 

Fj,{u)) as 



(J + J'-2)\ 



Fj,{lo) = 4 Cj_j, +2V /(2J+1)(2J' + 1) 



(J + J' + 3)! 
|(J||V 2 ||J')| 2 $>^ J '), (A10) 



v j' 



where n C 



~~ m\{n-m)\ 



represents combination. 



APPENDIX B: GEOMETRICAL FACTORS 



The geometrical factors Pf' for v — 3 and 4 in Eq. 
(|2.16|) have rather complicated forms. For v = 3, they 
are summarized as follows: 



expectation value of the rank-z/ multipole operator. 

J +2 



(A5) 



J' = J-2 



Ul - [-(2 J - 3)Fj_ 2 - (J — 2)Fj_! 

2Fj + ( J + 3)F/+i + (2 J + 5)F/ +2 ] , (A6) 
2" 



2 t /^[4(2J-3)(J-l)F 7 _ 2 
+(J-5)(J-l)f>_i 

-i(2J - 3)(2 J + 5)Fj + ( J + 2)( J + 6)F J+1 
+4(2J + 5)(J + 2)F J+2 ], (A7) 

2 [4(J-l)(2J-l)(2J-3)F/-2 

— (J — 1)(2J — 1)(J + 3)Fj_i 

+ (2J- 1)(2J + 3)Fj 

+ (J + 2)(J-2)(2J + 3)F J+ i 

-4(J + 2)(2J + 3)(2J + 5)F J+2 ], (A8) 

^[6J(J-l)(2J-l)(2J-3)Fj_ a 

-3J(J - 1)(J + 1)(2J - 

+J(J+ 1)(2J- 1)(2J + 3)Fj 

-3J(J + 1)(J + 2)(2J + 3)F/+i 

+6(J + 1)(J + 2)(2J + 3)(2J + 5)F J+2 ] . 



(A9) 



(3) _ : 



P 



p(3) 
A* 



1 



3\/2 



2 V 2 



p(3) 



[k x k] • q(e', e) + [e' x e] • q(k , k) 

+[k x e] • q(e', k) + [e' x k] • q(k, e) , (Bl) 
i [5 



(B2) 



[k x k^e^ + [e' x e^k'Jtp 
[k x el^kf, + [e x k] M fc^e M 
PjP for n = 2, 3, and 4, 



2V10 

j 13,,-/ 
4 



(4^ 



4 v 2 



7 

fl' ,fj," —'o 

(k'X6) M £ e W>"( e M'V - £ iL"V) 
7 



i /3 



7\/f( e ' x k) M £ f crt"(^e 



fx 1 ,fi"=5 

for = 5, 6, and 7. (B3) 



Note that /i = 2, 3, and 4 work as x, y, and z, respectively. 
Similarly, /1 = 5, 6, and 7 work as x, y, and z, respectively. 
The Levi-Civita tensor density e^'^" is introduced. 
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For v = 4, the results are as follows. 



P 



(4) 



P, 



(4) 



p(4) 



15 

p(0)i 



[5(k x k x e x e x + k k y e e y + k z k z e z e z ) 



(B4) 



p(4) _ 



2^^ [[k ■k] gi (e',e) + [e'-e]g 1 (k,k) 
+ [k • e]gi(e',k) + [e' • k] 9l (k , e) 

[k'-k] g2 (e',e) + [ e '-e]q 2 (k,k) 

+ [k • e]g 2 (e',k) + [e' • % 2 (k, e) 

1 6 
— j=Qn+i ( k . k ) XI e wV V - <V' 

1 6 

^— -=«Z^+i(e , e) ^ e w ' M » {k^k^i — k^, 

6 



(B5) 



1 



2V6 



(? M+ i(k,e) X e w v"( e M'^ 



M'./i"=4 
6 



(4) _ 



-i=g M+ i(e', k) X e w >»(^'V-^" 

2V0 M ', M »=4 

for /i = 4, 5, and 6, 

1 



(B6) 

m 

V') 
(B7) 



42 
1 



42 
1 



42 
1 



7fc p fc p 


-3(k 




_ 4 (e',e) 




- 3(e' • 




4 (k,k) 




-3(e' 


k)]g„- 


- 4 (k ,e) 




-3(k' 




-4(e',k) 



for ^ = 7,8, and 9, (B8) 



where indices 4, 5, and 6 in the summations in Eq. I|B7|) 
serve as x,y, and z, respectively. Similarly, indices 7, 8, 
and 9 appeared in the brakets in Eq. 1|B8|> serve as x, y, 
and z, respectively. 
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